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Summary—A survey of analytical models that predict the penetration of a shaped-charge jet into 
various target materials is presented. This survey discusses the early (1940 vintage) Bernoulli-like 
models, the extension of these models to include particulated jets and variable velocity jets and 
models dealing with compressibility effects and jet drift. Penetration models associated with metallic 
rods are also presented, especially insofar as they relate to shaped-charge or explosively formed 
penetrator studies. Modern analytical penetration models, which are used by the authors, are also 
discussed. An extensive bibliography is included. 


NOTATION 


Bmax  Brinell hardness of target 
C.D. charge diameter 
D rod diameter 
E, — energy in last part of projectile 
g gap distance between particles 
Jo empirical constant 
Hp dynamic hardness of target 
Hy height of interaction region 


k; jet body shape factor 

ky target body shape factor 
[jet length 

L rod length 

M _ jet mass 


P penetration depth 
Py penetration depth at zero standoff 
r, jet radius 
SorS.O. — standoff 
t time 
to time for jet to reach the target 
th jet breakup time 
T time at end of penetration 
U penetration velocity 
Unmin minimum penetration velocity 
V jet velocity 
V, velocity at the top of the interaction region 
Vinin Minimum jet velocity 
Vo jet tip velocity 
Y _ target strength 
a,B empirical constants 
y — (pr/p;)"? 
dg incremental length of jet 
A empirical constant 
~~ dynamic viscosity 
pj OF po jet density 
pr Or pt. _ target density 
a; dynamic strength of jet 
6p dynamic strength of penetrator 
Osp _ Stress state of deforming penetrator 
o, dynamic strength of target 
Other specialized symbols are defined where they are first used. 
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INTRODUCTION 


Analytical models capable of predicting the penetration of the jets from shaped charges into 
a variety of target materials are extremely valuable to the terminal ballistician. Many 
disciplines are involved, since shaped-charge jets are used to penetrate or perforate armor, as 
well as rocks, soil, wood, ice and other non-metallic materials. These analytical models 
provide fast, analytical predictions where the time and resources available prohibit large 
hydrocode computer solutions. Also, the target may be monolithic or it may consist of 
several layers of different materials including air or liquids. For these reasons, analytical 
penetration expressions usually assume one-dimensional flow and employ other simplifying 
assumptions. 

Early analytical penetration models (1940 vintage) were based on the Bernoulli principle. 
Later, empirical factors were introduced to account for particulation of the jet. Eventually, 
non-uniform-velocity (i.e. stretching) jets were incorporated into the models. Other authors 
included terms to account for jet and target strength effects, compressibility effects and the 
effects of jet drift, particle dispersion and particle tumbling. These models were based 
primarily on the Bernoulli principle or on a rod-type penetration model. Other models 
included transient effects based on a one-dimensional conservation of mass and momentum 
without recourse to the rod models or the Bernoulli concept. 


SHAPED-CHARGE JET PENETRATION 


Birkhoff et al. [1] and Hill et al. [2] developed a simple penetration theory from the 
hydrodynamic theory of impinging jets. This theory is documented in Refs [1-3]. 

Because of the hypervelocities associated with shaped charge jets, the pressures produced 
during jet—target impact far exceed the yield strength of most materials. Thus, to a first 
approximation, the strengths and viscosities of the jet and target materials can be neglected, 
allowing usage of the familiar hydrodynamic assumption of incompressible, inviscid fluid 
flow. 

Consider a shaped-charge jet of length /, density p; and velocity V penetrating a semi- 
infinite, monolithic target of density p;. The penetration velocity is U, as shown in Fig. 1. 
The penetration is simpler when viewed from a system of co-ordinates moving with the 
penetration velocity U, as shown in Fig. 2. In this system the hole profile is fixed, and the jet 
moves to the right at a velocity V — U and the target moves left at a velocity U. The pressure 
on the two sides of the interface between the jet and the target must be the same. Now, since 
the phenomenon is steady-state in this system of co-ordinates, Bernoulli’s equation may be 
applied along the axial streamline so that 


1 1 
5 e(V — U)? = 5 prU?. (1) 
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Fic. 1. Jet penetration. 
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Fic. 2. Jet penetration with a co-ordinate system moving at the penetration velocity, U. 


Figure 2 shows the jet being eroded while penetrating the target. Assuming that steady state 
is reached instantaneously and that the penetration ceases when the rear of the jet strikes the 
target, then the total penetration, P, is the penetration velocity times the penetration time or 
l 
V-U 


P=U or 


2 

P=I(p,/p;)'/?__ using equation (1). @) 

This simple equation indicates that the penetration is independent of the jet velocity. 

However, the rate at which the jet erodes and the final stretched length of jet depend on the 
jet velocity through the term /. Note that for real jets, / is a function of time. 

There are several limitations to this simple theory. Firstly, after the last jet particle strikes 
the target, i.e. after the jet is consumed, the residual inertia of the penetrator may cause the 
hole to continue to grow in depth (and width). This additional depth has been called the 
secondary penetration [3] or afterflow. The depth of the crater when the jet vanishes is called 
the primary penetration. 

Secondly, even the primary penetration is not predicted exactly, as effects other than 
density come into play. For example, a greater penetration into mild steel is observed than 
into armor steel, even though both have the same density. Moreover, more penetration into 
soft materials, such as lead, occurs than would be dictated by the density law. These results 
undoubtedly occur because material strength, strain, strain-rate dependence and other 
material properties are not included in the simple penetration model. These effects are 
especially important at low jet velocities, but less important at hypervelocities. In short, a 
density (and even a yield strength) is usually not sufficient to characterize a material. 

Thirdly, for real jets, the average penetration into a given target increases, reaches a peak 
value and then decreases as the distance between the base of the charge and the target (the 
standoff) increases. Figure 3 illustrates a typical penetration—standoff curve for a shaped 
charge with a conical liner penetrating armor steel. This variation with standoff is not 
predicted by the simple penetration model, as discussed later. 

If the jet is particulated (broken up) and separated such that the particles do not interfere 
with each other, the particles retain their velocity, V — U, and cross sectional area, A, until 
they impact the target, assuming the particles remain perfectly aligned and do not tumble, 
spread or diverge from the penetration axis. 

Upon impact, the dynamic pressure produced by a jet particle can be approximated by 
dividing the total force by the total cross-sectional area the jet strikes. 

Following Evans [4], the total force is given by the rate of change of momentum 
p;A(V — U)? and the average pressure on the impact surface is p(V — U)*. Equating this 
pressure to the pressure in the target material at the point of impact, from Bernoulli’s 
equation, yields 


1 
p{V—UY =; pyU?. (3) 


Comparison with equation (1) leads to 


Ap(V — U)? = p,U?, (4) 
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Fic. 3. Penetration—standoff curve for a conical liner. 


where 4 is a constant that equals one for a continuous jet and two for a particulated jet. p; is 
then the average jet density including the gaps between particles. For jets particulated 
partially or particulated during the penetration process, 1 < /4<2. Then from equation (2) 
the ideal primary penetration becomes 


1/2 
P= (2) c (5) 
PT 


This quantitative hydrodynamic theory of penetration by hypervelocity jets was 
developed by Hill, Mott, Pack and Evans (see [2], [4] and [5]) and independently by Pugh 
[6]. Additional details regarding the development of the hydrodynamic, hypervelocity 
penetration theory are given by Birkhoff et al. [1], Eichelberger [7], Birkhoff [8] and Ref. 
[3]. Birkhoff relates the parameter A to the drag coefficient and calculates the penetration for 
a rotating jet. 

This simple theory does not take into account jet velocity gradients, jet—target 
interactions, jet alignment, shock physics, compressibility effects, aerodynamic drag, 
variable-area jets, transient effects and jet-target material properties. References [5,7] 
provide additional information regarding penetration theory and secondary penetration (or 
afterflow effects). 

Other equations useful in penetration theory are, by definition, 


p= [var or U=P, (6) 
where P denotes the derivative of P with respect to time, and 
di =(V —U)dt. (7) 
Equations (6) and (7) yield 
U 
P= | V-u dl. (8) 


Early models attempted to improve the accuracy of the simple theory to better account for 
jet particulation (and spread) and standoff effects; these models were semi-empirical in 
nature. Let S denote the standoff and « be some constant depending on the jet velocity 
gradient; also, let 8 be a constant that determines the rate of jet spreading (dispersion). Then, 
for continuous (or very recently particulated) jets 


1+aS 
1+ BS’ 


where Py is the penetration at zero standoff distance [1]. For fully particulated jets where 
A =2, the penetration may be given as 


P= Pi[./2(1+a5S)"/2/(1 + BS)], (10) 


where P9 is the value of P from equation (9) at S = S,, where the jet breaks up [1]. This model 
assumes that the ratio of the effective jet length at some standoff § to that at S=Ois 1+4S. 


P=P, (9) 
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Also, the ratio of the effective cross-sectional area of the jet at standoff S to that at S =0 is 
assumed to be (1 + BS)? from [1]. Other models, due to Pugh [3, 6], assume that the ratio of 
the cross-sectional area of the jet at standoff S to that at S=0 is 1 + BS? instead of (1 + BS)?. 
In either model, the constants « and 8 must be determined experimentally. 

Pack and Evans [5] also noted the importance of target material strength on jet 
penetration. To account for strength, they proposed correcting equation (2) by a semi- 


empirical factor 
\1/2 Y 
p~(4) i(1- = i) (11) 
Pr pV 


They showed that for steel, the correction term aY /p,V? is as great as 0.3; that is, the effect of 
strength is to reduce penetration by as much as 30%. 

For the penetration of ductile targets, such as lead, by a shaped-charge jet, Pack and 
Evans [5] modified equation (11) by adding a secondary penetration term equal to the 
radius of the hole made by the jet. 

Significant progress was made in 1955, when Eichelberger [7] conducted penetration vs 
time measurements. He verified that the hydrodynamic formula, equation (1), is very 
accurate for early times and short standoffs where the jet has a high velocity and is not yet 
broken. He found that when the jet is broken, A in equation (4) should actually be less than 1 
if p; is taken as the original density of the liner. Furthermore, at low jet speeds, the strengths 
of the jet and target are important. Eichelberger proposed the formula 


Ap(V — U)? = p,U? +26, (12) 


where o = o, — g;,in which o, and g; are the ‘resistance to plastic deformation’ for the target 
and jet, respectively. These strength terms are taken as one to three times the static uniaxial 
yield stress, with the factor attributed to strain-rate effects and to the state of non-uniaxial 
stress at the yield site in the target. 

Beyond this point, the quasi-steady hydrodynamic theory has evolved along two lines: the 
shaped-charge jet with variable velocity, and the kinetic-energy rod with uniform initial 
velocity. For stretching jets, work has been directed toward analytical solutions and 
quantifying the degradation in penetration beyond the point of jet breakup. For rods of 
uniform velocity, models treating the deceleration and deformation (shortening) of the 
penetrator have been developed; these models are applicable to broken jets or short length- 
to-diameter ratio rods. 


Variable velocity jets 


For a jet of non-uniform velocity distribution, such as that produced by a shaped charge, 
the jet length is not constant but increases with time. In this case, equation (2) is not directly 
applicable. 

Asa first step, Abrahamson and Goodier [9] derived explicit formulae for the penetration 
of continuous straight jets of non-uniform velocity. They concluded that if strength is 
neglected, addition of a constant velocity over the entire jet actually decreases penetration, 
since the jet has less time to stretch before penetrating. If, however, the rear of the jet is too 
slow to overcome the strength of the target, addition of a uniform velocity may increase 
penetration. They derived a formula for the penetration of an idealized jet, which serves asa 
theoretical maximum. In their analysis, they took the jet length as known at a given distance 
from the target, but did not consider the relationship between jet length and distance from 
the charge; this makes their results less useful for practical purposes. 

Allison and Vitali [10] further extended the theory to account for jet segmentation, with 
the following assumptions and conclusions. 

(a) Existence of a virtual origin. It is known that the velocity of each jet particle remains 
nearly constant. If, in addition, the spatial distribution of velocity is linear, as for most jets 
from constant-thickness conical liners, then a ‘virtual origin’, or point source, of the entire jet 
can be located, as first proposed by Allison and Bryan [11]. 
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(b) Negligible strength. The strengths of the jet and target materials were neglected by 
Allison and Vitali [10]. So, to characterize the termination of the penetration by the slow- 
moving rear of the jet, a minimum jet velocity for penetration, V,,,,, is defined. This V,,,,, and 
the penetration process at low velocity, must depend on the strengths of the jet and target. 
This is an artificial method to terminate the penetration. 

(c) Negligible compressibility. They also studied the effect of compressibility on 
penetration from equation-of-state data and the compressible-flow Bernoulli’s equation. 
They concluded that firstly there is no effect for like jet and target materials and secondly, for 
dissimilar materials, the effect is small, provided the compressibilities do not differ greatly 
(e.g. only 1°% difference in penetration velocity between the compressible and incompressible 
flow equations for a copper jet against a titanium target at shaped-charge velocities). Harlow 
and Pracht [12] confirmed this in two-dimensional code simulations of iron and aluminum 
jets impacting similar targets. (It will be shown below that compressibility is indeed 
negligible for a metal jet on a metal target; for a metal jet on a low-sound-speed target, such 
as Plexiglas, compressibility is not insignificant.) 

Before jet breakup, penetration P as a function of time t was derived by Allison and Vitali 
as 


P(t) = Vot(to/ty/"*” — toVo, (13) 


where » = (p,/p;)'/? and fy is the time of arrival at the target of the jet tip, which has velocity 
Ve. 
(d) Simultaneous breakup. In Ref. [10], the entire jet is assumed to break simultaneously. 
This assumption has also been used by other investigators [13-17]. In general, for 
conventional shaped charges, where the jet particulates or breaks up from tip to tail, the 
simultaneous jet breakup assumption is, at best, a first-order approximation. Usually, as 
verified by the studies of Simon et al. [13-17], this assumption is satisfactory for 
conventional warheads. 

In reality, there exists a distribution of jet particulation from tip to tail where several 
microseconds are required before particulation is complete. For unconventional warheads, 
where the breakup does not proceed from tip to tail, this assumption (of a simultaneous, 
constant, breakup time) is invalid. 

(e) Each broken-jet segment penetrates as a continuous jet. The basic principles involved in 
Eichelberger’s hydrodynamic penetration theory are Bernoulli’s equation and the continuity 
equation, 


Se (14) 


where dé is the incremental length of the jet. Allison and Vitali, in calculating the penetration 
of a segmented jet, assumed that equations (1) and (14) are applicable for short separated jet 
segments as well. This assumption is questionable, because the penetration of finite segments 
is affected more by unsteady effects and is often less than that calculated from steady 
formulae. Under this assumption, once the jet breaks, penetration remains constant with 
further increase in standoff. In reality, penetration decreases with standoff after the jet is 
broken. Experiments by Chou et al. [18, 19] indicated that the penetration of two segments 
with an intervening gap is less than that of two segments with no gap. Two-dimensional 
computer-code calculations by Harlow and Pracht [12] also demonstrated that, during the 
initial transient period, penetration is less than in the steady-state case. 

The principles proposed by Allison and Vitali are still used today as the basis for 
penetration calculations. Dipersio and Simon [14] presented explicit formulae based on 
Allison and Vitali’s theory, for three cases: (a) penetration before jet breakup (T < t,), (b) jet 
breaks during penetration (t9<t,<T) and (c) jet breaks before reaching the target 
(t, <to < T), where T is the time at the end of penetration, f, the time of jet breakup and ft) 
the time when the tip reaches the target. For case (a), the total penetration depth P is 


P =S[(Vo/Vein)'"” — 1], (15) 
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where S here is the distance from the virtual origin to the target, or the effective standoff 
distance. In this case S is bounded by 


V...\t/ 
0<S< Vas ss) 
Vo 
Vein 


l+y 


and ae = 


For case (b), the penetration is 


palit rity) Mt #56 — Vent 


S, (16) 
y 
where S is bounded by 
Vinin ae 
Fonts V, ) <S<Vot, 
Finally, for case (c), the penetration is 
Vo — Vouin)t 
P= ( 0 min) b (17) 
? 
for standoffs in the range given by 
Votp<S<m. 


Recall that y = ./p7/p;. 


The formulae given above are the results of the so-called DSM (DiPersio, Simon and 
Merendino) theory [14]. Other works related to studies of this type are given in Refs [13-17, 
20-26]. 

The above equation may be re-arranged to yield the velocity V, of that point in the jet that 
is currently penetrating at a depth P in the target. For example, equation (15) becomes 


S y 
= Hl prs | | 


This is also, of course, the emergent velocity of a continuous jet after penetrating a target of 
finite thickness P. 

Based on equations (15-17), DiPersio and Simon [14] also compared their calculated 
penetration depth with experimental data. In general, good agreement was obtained at short 
standoffs, up to about three charge diameters. At larger standoffs, the formulae give higher 
than measured values. They attributed the declining performance of actual jets mainly to 
asymmetric wavering of the jet. 

In 1965, DiPersio, Simon and Merendino [17] conducted an extensive experimental 
program to determine the causes of decreased penetration at long standoff. They reached 
the following conclusions: 

(a) Jets from precision charges travel in very straight lines; jet particles show no 
appreciable tumbling, deceleration or mass ablation due to air friction. For conical liners, 
the jet-velocity distribution is nearly linear, and the assumption of a virtual origin is valid. 

(b) The minimum jet velocity for penetration is not constant for a given jet or target, but 
increases with standoff. 

Since V,,;, is not a constant, they proposed a new criterion for terminating penetration 
based on a minimum penetration velocity U,,,. From test data they noticed that the 
penetration velocity U at the time penetration stops was constant for different standoffs. The 
penetration velocity of a broken jet was obtained, in essence, by drawing a smooth curve 
through the penetration history. The slope of this curve at the cessation of penetration was 
taken to be a constant, U,,,,, regardless of standoff. For the standard charge they used, U,,., 
was about 1.0kms™!. 
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Next, DiPersio and Simon modified the penetration formulae in all three regimes by 
replacing the V,,;, criterion with a U,,;, criterion. The new formulae yield penetrations that 
decrease with standoff, in better agreement with experimental data [17]. 

Equations (15-17) were taken to be valid except those equations that contain the term 
V.jn- Equations for the jet velocity and the penetration velocity were derived. The equations 
for the total penetration, using U,,;, aS a jet cut-off criterion instead of V,,;,, are given below. 


For case (a), 
7 Vo ly 
P=S(azeuc) 7 | ne 


where S is bounded by 


1+ y)Umig ]!” 
0<S< (+ )Umah| ODE e | ‘ 
0 


For case (b), the penetration is now 


pt (Vor) PSHE? — J+ Wrsintel Voto)! 8 *S0*? 


vw 
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S, (19) 


where S ranges from 
(1+ y)U nin 


1/y 
<S< Voty. 
0 


(1 + Want 


For case (c), when the jet particulates before penetration begins, 


P= Voty — ¥/ ae + 7S) (20) 


y 
where S lies in the interval 


a 


) 
f 


Votp [Vo — Umi 
Votp<S< fs (To Sa) 


min 
and P = 0 when S equals its upper bound. 

For the U,,,, criterion, equations (18-20) are used in lieu of equations (15-17), 
respectively. The total length of jet which contributes to the penetration process may be 
obtained by multiplying the appropriate value of P [from equation (18), (19) or (20)] by y. 

DiPersio and Simon subsequently conducted many informative experiments. In 1966 
[15], they reported a study on ‘residual penetration’—penetration by a jet into a main target 
after perforating a series of skirting plates. In 1968 [16], they studied the effect of jet breakup 
on penetration (the earlier the breakup, the less the penetration), and discussed the effect of 
jet rotation. In Refs [20, 23, 26], they verified experimentally that penetration decreases as 
target strength, or hardness, increases. Merendino and Vitali made the same observation 
[21]. 

In Ref. [22], Simon and DiPersio, using flash radiography and other techniques, revealed 
that the rear portion of the jet that does not contribute to penetration actually piles up at the 
bottom of the hole. They also mentioned a ‘valve action’, whereby the rear of the jet is 
trapped by the small hole created by the front of the jet. References [17, 20, 22-26] present 
several comparisons of the DSM theory with experimental data. 

Note that in all these studies the basic hydrodynamic equation (1) was used to describe the 
penetration process, while efforts were made to determine the termination of penetration, 
Vain OF U,nin, by empirical means. 

Eichelberger’s formula, equation (12), which contains the material strength term, was not 
used at all. Thus, the material is described only by an empirical constant and its density. 

The penetration criteria based on the minimum velocities V,,,, and U,,,, provide useful 
practical formulae, but do not directly address the physical causes of the decrease in jet 
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penetration with standoff. These phenomena, which include jet segmentation, wavering and 
tumbling of jet particles, have been the focus of several recent studies. Note that since U is the 
penetration velocity and by definition U = P, so U,,;, (the minimum penetration velocity) is 
zero when the penetration ceases, or when the penetration—time curve asymptotes. 

It is generally recognized that three-dimensional effects (or two-dimensional effects if the 
target is impacted at normal obliquity) are the true mechanisms involved in jet cutoff or 
penetration termination. A U,,,, OF Vi, cutoff criterion is simply a device to allow the one- 
dimensional models to be effective. This is not meant asa criticism of the DSM theory alone, 
since all analytical models suffer this limitation. Basically, the DSM limitation is in the a 
priori determination of U,,;,, which is often not a constant, but varies strongly with standoff 
and charge diameter for a given jet—target configuration. This is documented in the CUMIN 
report [27] where Majerus and Scott provided several plots of the U,,,, variation. 
Unfortunately, this variation of U,,;, with standoff and charge diameter is highly erratic and 
cannot be readily modelled or curve fit. Nonetheless, Ref. [27] is a useful aid in the 
application of the DSM theory. 

Jet wavering, jet tumbling, jet drift and crossing (or transverse) velocity effects (where the 
target moves as the jet penetrates) have been analyzed by several authors. Notably, these 
analyses are found in Majerus et al. [28] and Segletes [29] on the effects of transverse 
velocity, Segletes [30] on the effects of jet drift velocity and Held [31] and Golesworthy [32] 
on the effects of transverse velocity. 


Particulated jets 


Carleone et al. [33], in an augmented version of the one-dimensional jet-formation code 
DESC, developed a theory for predicting the breakup of shaped-charge jets and introduced 
an approximate method for determining the decrease in penetration with standoff. First, the 
breakup time of each portion of the jet is calculated or measured. If a given segment of jet 
reaches the target before its computed breakup time, its penetration is calculated from 
Eichelberger’s equation (4), for a continuous jet. If the segment of jet has broken, however, its 
actual penetration dP’ is deduced from the continuous-jet penetration dP by the formula 


ap’=<r{ 1 ~£), (21) 
Go 


where g is the gap distance between jet particles, non-dimensionalized with respect to the 
increment of jet length, and gy is an empirical constant. According to this formula, 
penetration is degraded as the particles move farther apart, so that after breakup, 
penetration decreases (as expected) with standoff. Comparisons with experimental data 
show good agreement when gy = 6.5 for precision shaped charges and 4-6 for non-precision 
or small charges. An advantage of this simple theory is that only one empirical constant need 
be specified. 

A more sophisticated approach was taken by Smith [34], who applied the Monte Carlo 
technique to the effects of tumbling, dispersion (or wavering) and breakup, which are 
considered as stochastic processes. The values of the parameters governing these processes 
are assumed to follow prescribed probability distributions. Then, random values consistent 
with these distributions are chosen for the parameters of each jet particle. Finally, jet 
penetration is determined by applying Tate’s method for rod penetrators [35] (discussed 
later) to each particle. With proper selection of the statistical parameters (mean, standard 
deviation) describing the distributions of tumbling rate, dispersion angle and breakup time, 
excellent agreement with measured penetrations can be achieved. The chief shortcoming of 
this approach is the difficulty in determining, a priori, proper values for all of these 
parameters. Nevertheless, this method is useful for interpolating and extrapolating 
penetration data; once the parameters have been calibrated by experimental data, 
penetrations at other standoffs or against other targets can be estimated with some 
confidence. Improvements to this model [34] are possible by coupling the Monte Carlo 
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technique with a model more suitable to shaped-charge penetration than Tate’s rod model 
[35]. 

Fa [36] formulated a penetration—perforation model based on the conservation of mass 
and momentum. Zaid considered shear stress terms and modelled a non-deforming 
cylindrical penetrator. 

Majerus and Walters [37-43] developed a penetration model based on the conservation 
of mass and momentum using an integral approach and a control-volume concept. They 
divided the penetration process into two regions, the undeformed portion of the penetrator 
and the interaction region of deforming target and penetrator materials at the bottom of the 
crater. Global equations of motion are derived for each region in terms of normal stresses 
and shear stresses, the latter using a Newtonian fluid approximation. The resulting system of 
ordinary differential equations is solved by numerical integration. With proper selection of 
values for target strength and viscosity, good agreement with measured penetration can be 
achieved. The dynamic viscosity values are obtained from the U.S. or U.S.S.R. experimental 
data. See Refs [36-42], for example. The strength term represents a dynamic value 
analogous to that used in equations (11) and (12), for example. 

Walters and Majerus used an approximation of broken jets involving two steps. First, the 
broken jet is replaced by a continuous jet of equal length, of density the same as the average 
of the broken jet, as suggested by Allison and Vitali [10]. Then the continuous jet is 
shortened to the sum of the lengths of the particles in the original jet. This method yields 
satisfactory results but gives no insight into the mechanism of penetration by broken jets. 

The model defines the termination of penetration to occur when U, the penetration 
velocity, reaches zero. This can happen after the jet length (which is also calculated) is 
consumed, thereby providing a measure of the secondary penetration or afterflow. The 
model is also capable of analyzing layered, multi-material targets including air gaps [41]. 

The conservation of momentum yields 


“ (MU) = ar?(o, — 07) + ar? p9(V; — U) — ar}? pz,U? — 2nH wU (22) 


for a continuous jet. In this equation M is the mass of the jet, U is the penetration velocity, 7; 
is the jet radius, o, and o; are the dynamic strength values of the jet and target, respectively, 
Po and p;, are the jet and target densities, respectively, Ho is the height of the interaction 
region between the penetrator and the target and V, is the velocity at the top of the 
interaction region. References [37—43], especially [41], present the derivation of equation 
(22) and provide analogous equations for particulated jets. Expressions for V,, Ho, r, and the 
jet length are also given for continuous and particulated jets. The jet length, jet radius and jet 
density determine the jet mass M. References [37—43] present several penetration depth and 
penetration time predictions and comparisons with experimental data. References [37-40] 
extend the model to include penetration by long or short rods. 

One interesting application of this model was the accurate prediction of the data presented 
in Ref. [3]. In that work, a series of 1 in. plates of lead were stacked and a single | in. 
RHA (rolled homogeneous armor) plate was inserted in the stack at various depths, e.g. first 
as the top plate, then as the second plate, etc. The test objective was to plot penetration into 
the stack as a function of the position of the hard (RHA) plate from the top of the stack. For 
penetration by a shaped-charge jet, there exists an optimum position of the hard plate from 
the top to minimize the penetration. The experiment was also conducted using mild steel 
plates instead of lead. The Walters-Majerus model [43] successfully predicted the optimum 
location of the RHA plate for either case given in Ref. [3] and for a new experimental test 
using an RHA plate in an aluminum stack. In addition, this effect was found to be standoff- 
dependent; see Fig. 4, for example, where an RHA plate in a lead stack was analyzed. This 
study was significant in that many models experience difficulty in providing accurate 
penetration depth and penetration time predictions into multi-material layered targets since 
any material parameters or semi-empirical constants in use must change during the 
penetration process. 
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Fic. 4. Predicted penetration vs a layered target, from Walters et al. [43]. 


Most other models in use are based on an extension or duplication of those described in 
this section. For example, a simple analytical model based on the concepts given earlier was 
developed by Szendrei [44]. Other models that contain original concepts are due to Evans 
and Ubbelohde [45, 46], Luttwak et al. [47], Andersson et al. [48], Sagomonyan [49] and 
Murphy [50]. 

A frequently used model for shaped-charge jet penetration follows from Alekseevskii [51], 
Sanasaryan [52] and Sagomonyan [53-55]. The simple penetration model is based on a 
modified Bernoulli equation, 


Hp + ky pp U? = asp + k,p,(V — U)’, (23) 


using the notation introduced earlier, and where Hp is the dynamic hardness of the target, 
taken as the Vickers hardness. The terms k; and k; are the body shape factors of the target 
and jet, respectively, which depend on the flow geometry and the effect of internal forces; 
they are difficult to determine—Alekseevskii [51] takes k; = k; = 1/2. Finally, asp is the part 
of the stress state of the deforming penetrator determined by the internal forces, and is called 
the dynamic yield point [51]. The cutoff velocity, or velocity below which the jet will cease to 
penetrate, akin to the V,,,, of the DSM theory discussed earlier, is 


Vern = [Ap ~ %sp. (24) 
kp; 


Equation (24) follows from equation (23) for U = 0. Since U = P, determination of U implies 
P by numerical integration. 


Compressible models 


As discussed earlier, Allison and Vitali [10] and Harlow and Pracht [12] indicated that 
the effect of compressibility on penetration of metal jets into metal targets is slight. However, 
Haugstad and Dullum [56, 57] showed that for non-metallic targets, compressibility effects 
may be significant. Furthermore, if the jet velocity is sufficiently high for shocks to occur, the 
penetration can be much less than that predicted by the incompressible theory (20% 
reduction for a copper jet against plexiglas at 10kms~?). 

The form of the Bernoulli equation presented earlier neglects compressibility, 
consideration of which requires an additional term (the flow work) that depends on the 
compressibilities of the jet and target materials. If the compressibilities are similar (as for a 
copper jet impacting a steel target), this term is negligible. If the target is appreciably more 
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Fic. 5. Effect of compressibility on penetration by a uniform-velocity copper jet [58]. 


compressible than the jet (as for a copper jet against a plexiglass target), the effect of this term 
is to significantly reduce penetration as compared to incompressible flow theory, as shown in 
Fig. 5. 

Flis and Chou [58] extended the compressible theory to include more general equations of 
state (including Mie-Gruneisen and Tillotson), and reconfirmed the conclusions of previous 
researchers. They also applied the compressible model to a stretching jet using the virtual- 
origin assumption, and showed that the penetration of such a jet can theoretically be as 
much as 50% less than predicted by the incompressible theory. This greater degradation in 
penetration for a stretching jet is explained by observing that, when the front of the jet 
penetrates less, the rear portion has less distance to stretch, and therefore also penetrates less. 

Experimental evidence of the importance of compressibility has not been convincing. 
Experiments by White et al. [59,60] of non-precision shaped charges against plexiglas, 
aluminum, steel and several liquids as target materials showed no correlation of target 
penetration resistance with compressibility. 

Chick et al. [61] applied the compressible theory to a jet with non-linear velocity 
distribution, for which a virtual origin was not assumed. For such a jet projected by a small 
copper-lined conical charge into a plexiglas target, a difference in penetration rate of only 
about 4% is predicted by the compressible and incompressible theories, but they claimed 
that their experimental results were in better agreement with the compressible theory. 
Coincidentally, they point out that, for their jet, the error involved in neglecting 
compressibility is comparable to that involved in assuming a single virtual origin for the 
entire jet. 

Woidneck [62] fired copper rod projectiles into water, glycerol and methanol at velocities 
of 2500 and 5000 m s~'. Even at the higher velocity, the effect of compressibility makes a 
difference of only 3.3% (for water) in the early-time penetration rate, which he did not 
discern in his experiments. 


The virtual origin concept 
The virtual origin concept is used in penetration calculations to determine the effective 


standoff distance and to allow an accurate prediction of short standoff penetration. Most of 
the penetration models assume that the jet velocity is linearly distributed in space and rely on 
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the virtual origin approach first used by Allison and Vitali [10] to integrate the penetration 
over the length of the jet. Recently, Chou et al. [63] examined the accuracy of the virtual 
origin approximation, concluding that for charges with simple conical or hemispherical 
liners, the maximum error involved is less than 0.2 charge diameters and is constant with 
standoff. However, Foster [64], using a mass-summation technique to analyze jet 
radiographs before breakup, has shown a significant non-linearity in the velocity 
distribution of a biconic shaped charge. 


Uniform velocity rods 


Uniform velocity rod models are often used in shaped-charge penetration and SFF or 
EFP penetration studies. Note that there is some analogy between penetration by rods and 
shaped-charge jets in that pure rod models are sometimes adapted for shaped charge 
penetration studies, e.g. Smith [34], Alekseevskii [51] and Sanasaryan [53]. Some 
noteworthy differences are in the target strength effects which are usually more predominant 
for rod-like projectiles which experience typically lower impact velocities than shaped- 
charge jets. Thus, penetration by a shaped-charge jet represents a closer approximation to 
purely hydrodynamic flow than rod projectiles. For a shaped-charge jet, the jet velocity can 
be established, since the velocity gradient of the jet from tip to tail is approximately linear. 
However, for rod-like projectiles, which deform and erode and do not penetrate at a constant 
velocity, the rod velocity during penetration (as well as the penetration velocity) is unknown. 
On the other hand, for shaped-charge jet penetration studies, one must consider the eventual 
breakup of the jet and analyze penetration of segmented jet particles. 

Allen and Rogers [65] conducted a series of experiments with rods of various materials 
impacting aluminum targets. Comparing their data with Eichelberger’s formula, equation 
(12), they concluded that the net strength ¢ = o, — a, is nearly independent of penetrator 
material, or that o, is negligible compared to o,. Further, they postulated that the target 
strength is strain-rate dependent and expressed o, as a function of the penetration velocity U. 
In addition, they observed ‘secondary penetration’ in high-velocity impacts of dense (gold) 
rods against light (aluminum) targets. In this case, at the end of the usual penetration 
process, penetrator material that has been deflected by the crater bottom still has a net 
velocity into the target and can thus cause additional penetration. 

Christman and Gehring [66] developed a model by considering the effects of the primary 
and secondary phases of penetration. For a velocity range of 2.0-6.7 kms~', they report 
good agreement with experimental data. The dimensionless penetration into semi-infinite 
metal targets is expressed as 


P 1/2 D 2/3/9 y2\1/3 
Pa(1—-2\(Pe) 42.42? (Pe) (Ar) (25) 
L L Pr L Pr Braax 


The first term represents the primary (or hydrodynamic) phase of penetration. The effective 
length of the penetrator for this phase is reduced by one rod diameter. This remaining 
diameter of length is assumed to contribute to the secondary phase of penetration, as 
represented by the remaining empirical term. The terms are self-explanatory, except perhaps 
Byax, the Brinell hardness of the target. This model allows the final penetration to be greater 
than that predicted by the one-dimensional, hydrodynamic theory. 

Doyle and Buchholz [67] used a modified form of the Christman—Gehring equation (25) 
to predict the penetration of mass focused slugs (or EFPs) at long standoff distances. Their 


equation was 
P _(1_D\(ee)"” 043 (o0)""( Bs \"” 06) 
L L Pr L Pr Brnax : 


P = target penetration depth (inches), 
L=length of projectile (inches), 
D =diameter of projectile (inches), 


where 
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pp = density of projectile, 

p, = density of target, 

E, = energy in last caliber (last part) of the projectile (joules), 

Buax = target hardness beneath the penetrating projectile (kg mm 7’). 

From Ref. [67], Doyle and Buchholz state that equation (26) essentially assumes that the 
first part of the slug (all except the last caliber) performs like a shaped-charge jet. The last 
part (last caliber) performs like a ballistic projectile. Thus, total penetration in metallic 
targets can be represented by this two-part equation that separates the contributions from 
primary and secondary penetration phases. Primary penetration is a function of projectile 
length and density and target density. Secondary penetration is a function of projectile 
density and kinetic energy and target density and strength. 

A major contribution to the one-dimensional theory of rod penetration was made by Tate 
[35, 68-73], who accounted for the deceleration and shortening of the rod during 
penetration. (A similar model was developed independently by Alekseevskii [52].) Equation 
(12), the Bernoulli equation with strength, is used with the equation of motion of the 
undeformed part of the rod 


(27) 


and the kinematic relation dL/dt = V — U, where Lis the current length of the rod. Tate 
presented several explicit solutions to this system of equations for special cases of density and 
material strength; the general case may be solved by numerical integration. Tate 
recommended the value of c, as the Hugoniot elastic limit of the projectile material and o, as 
3.5 times the Hugoniot elastic limit of the target material. 

Tate considered two classes of impacts, o, > 0, and o, <a,. In the former case (hard 
penetrator, soft target), the penetration process passes through two stages: at a sufficiently 
high impact velocity, the penetrator and target first both flow hydrodynamically; later, after 
the penetrator has decelerated, it ceases to flow and penetrates as a rigid body. When a soft 
penetrator impacts a hard target (the case of most interest), both materials initially behave 
hydrodynamically; then, after the penetrator decelerates below a critical velocity, the target 
ceases to flow so that penetration stops, while the penetrator continues to deform. For both 
cases, if the impact velocity does not exceed a critical velocity calculated from equation (12), 
the initial stage of purely hydrodynamic behavior is absent. These phenomena are 
summarized in the phase diagram, Fig. 6, introduced by Chou and Flis [74]. 

Chou and Flis suggested a different set of values for o, and o,. They reasoned that, since 
the deforming part of the penetrator is in a state of uniaxial stress, ¢, should be chosen as the 
uniaxial yield stress rather than the Hugoniot elastic limit (the strength under a condition of 
uniaxial strain). Then, by comparing with experiments and hydrocode calculation, they 
showed that taking the target resistance o, as 2.5 times the Hugoniot elastic limit gave the 
best agreement. In addition, in a computer code for integrating Tate’s equations [74], they 
included an approximate method for computing the secondary penetration observed by 
Allen and Rogers, and also by Tate. 

An Integral Theory of impact has been developed by Donaldson and Contiliano and 
others [75,76]. This theory was extended to include a rod penetrator model applicable to a 
wide range of impact velocities and materials. Good correlations with experimental data 
have been reported. The penetrator is modelled as a head, from which material can erode, 
and a rigid shaft that is decelerated by interface stresses at the head. The target material is 
characterized by two properties: E,,, the energy absorbed during plastic deformation, and 
E,,.., that absorbed by elastic deformation. Simple formulae have been derived to determine 
these constants from laboratory tests. The determination of these constants for a variety of 
materials has identified certain hard, brittle, non-metallic materials as effective armor 
materials. 

In Fig. 7, three calculations are compared based on Tate’s deceleration model (using the 
RODPEN code [74]), Donaldson’s integral theory of impact and Christian and Gehring’s 
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Fic. 6. Phase diagram of rod penetration according to one-dimensional theory. 
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Fic. 7. Comparison of three one-dimensional theoretical calculations with experimental 
penetrations of steel rods impacting armor steel. 


semi-empirical formula with impact experiments performed by Hohler and Stilp [77]. These 
results show good agreement among RODPEN, Donaldson and the experiments, while 
Christman and Gehring’s formula seems useful here only at the higher velocities, and indeed 
was developed from high-velocity data. 

Matuska [78] investigated steady-state jet penetration using a computational approach. 
From simulations with the two-dimensional Eulerian hydrocode HULL, he determined the 
values of the parameters in a modified form of Bernoulli’s equation, 


p,U? 
3 


5 AV ~U)? + Bo, =~ + a0, (28) 
The values determined for the parameters are: 

a=1 

B=0.3 

y = 0.47 + 0.0289; + 0.00086p? + 0.072 In V. 
The deviation of y from its theoretical value of unity was attributed to the target’s resistance 
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to radial flow. This deduction was based on results of simulations of two impinging jets of 
equal radius, for which » had a value very close to 1. 

This model was extended to rod penetrators by substituting the above equation for the 
Bernoulli equation in Tate’s model. The steady-state portion of penetration is calculated by 
numerically integrating the resulting equations. Steady-state penetration is considered to 
cease when the rod has eroded to a length equal to its diameter (as suggested by [66])} or 
when the penetration velocity has dropped to a prescribed value. 

Other rod models and data are discussed by Walters and Majerus [37], Belyayev et al. 
[79], Weihrauch [80] (an excellent experimental study), Clark [81] and Dehn [82], who 
discusses the history of rod penetration, presents a survey of rod penetration models and 
develops new rod penetration concepts based on estimated or semi-empirically determined 
parameters. Allen et al. [83] introduce a drag coefficient concept and present a model for a 
non-eroding bullet penetrating sand. 

Excellent surveys of penetration by rods are given by Wright [84], Backman and 
Goldsmith [85], Dehn [82], and in the books by Goldsmith [86] and Johnson [87]. Studies 
involving rod penetration, experimental data and hydrocode simulation are given by Jonas 
and Zukas [88,89] and Zukas [90]. Zukas summarizes much of the penetration and 
hydrocode simulation methods in his book [91]. 

As a final note, Chou and Flis [92] recently published, in part, some of the material 
presented earlier, such as one-dimensional hydrodynamic theory, two-dimensional theory, 
variable-velocity jets, compressible penetration models, virtual origin concepts and uniform- 
velocity rod penetration. The current report elaborates upon and expands the abbreviated 
discussion of [92]. 

Also, Dehn [93] recently published a generalized equation of motion for penetration. In 
various specialized forms, this equation is used to calculate penetration by rods, jets or 
fragments. The generalized equation of motion describes the target force resisting 
penetration as a second-order polynomial in the penetration speed. The force and 
coefficients of the penetrator speed must be determined for each case. This equation was also 
presented by Johnson [87] and Allen et al. [83] and is a generalized form of the 
Robins—Euler equation. Dehn [93] also provides a comprehensive overview and survey 
of jet and rod penetration models. 


CONCLUSIONS 


Several analytical models used to calculate the penetration of a shaped-charge jet were 
presented. These models range from those of historical value, which provided the foundation 
for later models, to models currently in use. An overview of some of the available models was 
presented. No comparison was made regarding the relative accuracy of any of these models 
since they all involve either dynamic material properties of unknown value or semi-empirical 
parameters that must be determined for the given penetrator—target configuration. Thus, the 
choice of a penetration model depends on the input available to the user and his experience 
with a particular model. 

An extensive reference list is provided to guide others in choosing a model or in developing 
a new model. 
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